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Abstract 

We quantitatively explore, in the context of the DO detector at the Tevatron, three 
very different techniques for observing delayed decays of the lightest neutralino of a simple 
gauge-mediated supersymmetry breaking (GMSB) model to photon plus gravitino. It is 
demonstrated that the delayed-decay signals considered can greatly increase the region 
of general GMSB parameter space for which supersymmetry can be detected. In the 
simple class of models considered, the combination of standard supersymmetry signals 
and delayed-decay signals potentially yields at least one viable signal for nearly all of 
the theoretically favored parameter space. The importance, for delayed-decay signal 
detection, of particular detector features and of building a simple photon detector on the 
roof of the DO detector building is studied. 



1 Introduction 

Even in the context of the Minimal Supersymmetric Model (MSSM), it has become increas- 
ingly apparent that there are many different ways in which supersymmetry breaking can occur, 
leading to many different possible phenomenologies. In particular, gauge-mediated supersym- 
metry breaking (GMSB) models have attracted a great deal of interest in the last few years. 
The theory of gauge-mediated supersymmetry breaking posits that supersymmetry breaking 
is transmitted to the supersymmetric partners of the Standard Model (SM) particles via the 
SU(3)xSU(2)xU(l) gauge forces. Two GMSB model-building approaches have been explored 
in the literature. 

• In hidden-sector models, the GMSB model consists of three distinct sectors: a hidden 
sector where supersymmetry is broken, a "messenger" sector containing a singlet field and 
messenger fields with SU(3)xSU(2)xU(l) quantum numbers, and a sector containing the 
fields of the MSSM The coupling of the messengers to the hidden sector generates 
supersymmetry-breaking in the messenger sector. 



• In models of "direct gauge mediation" (sometimes called direct-transmission models) 
0, the GMSB model consists only of two distinct sectors: a sector which is not only 
responsible for supersymmetry breaking but also contains the messenger fields; and the 
sector of MSSM fields. 

In both classes of models, supersymmetry-breaking is transmitted to the MSSM sector via the 
SM SU(3)xSU(2)xU(l) gauge interactions between messenger fields and the MSSM fields. In 
particular, soft-supersymmetry-breaking masses for the gauginos and squared-masses for the 
squarks and sleptons arise, respectively, from one-loop and two-loop diagrams involving the 
virtual exchange of messenger fields. 

Let us consider further the simple hidden-sector models of the first class in which the 
communication of supersymmetry breaking in the hidden sector to the messenger sector occurs 
via two-loop interactions involving a new gauge group with coupling Qm and in which the 
messenger sector consists of an SU(5) singlet superfield S and a certain effective number Nm 
of complete SU(5) representations. In such models the boundary conditions for the soft- 
supersymmetry-breaking masses are set at the mass scale of the messenger sector Mm and 
take a particularly simple form. For the gaugino Majorana mass terms one finds [|T|: 

M,(M^) = hNmg (^) , (1) 

where ^2 = ^3 = 1,^1 = 5/3 and Nm < 4 is required to avoid Landau poles. A = {Fs) / (S) 
is the ratio of the vacuum expectation values of the auxiliary and scalar field components of 
the chiral field 5*. Note that the gaugino masses have the same relative values as if they were 
unified at the grand-unification scale Mu even though the actual initial conditions are set at 
scale Mm- In addition, scalar masses are expected to be flavor independent since the mediating 
gauge forces are flavor-blind. In the simple models we consider here, they take the form: 

^nM^).2A^N„j{A.)p..[<^J (2) 

with C3 = 4/3 (for color triplets), C2 = 3/4 (for weak doublets) and ci = | (^^^ (in the 
normalization where Y/2 = Q — T3). To avoid negative mass-squared for bosonic members of 
the messenger sector Mm /A > 1 is required; Mm/ A > 1.1 is preferred to avoid fine-tuning, 
for which /(A/M^) ~ 1 and 1 < g{A/Mm) < 1.23. For Mm/A > 2, 1 < g{A/Mm) < 1.045. 
Here, we will consider masses obtained with g = f = 1. The results above at scale Mm must 
be evolved down to the scale of the actual sparticle masses, denoted Q. The resulting gaugino 
masses are given by replacing ai{Mm) in Eq. (|I]) by ai{Q). Evolution of the sfermion masses 
is detailed in Ref. |]T^. Most important is the ratio mj^/Mi = \jQ/^\Jri/Nm — (5/33)(l — ri), 
where ri = {ai{Mm) / O-iiQ))"^ ■ For the very broad range of Q > mz and Mm < 3 x 10^ TeV, 
1 < ^1 ^ 1-5, in which case the lightest of the sparticle partners of the SM particles is the 
x\ ^ B for Nm = 1 or the in (more precisely, the fi) for Nm > 2. Our numerical results are 
obtained by evolving from Mm = 1.1 A. 
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As for the other parameters of low-energy supersymmetry, the A-parameters (describing 
tri-hnear Higgs-squark-squark and Higgs-slepton-slepton couphngs) are suppressed, while the 
generation of the /i and B parameters (responsible for mixing of the Higgs superfields and scalar 
fields, respectively) is quite model-dependent (and lies somewhat outside the standard ansatz of 
gauge- mediated supersymmetry breaking). In the end, the parameters of the minimal GMSB 
model P] consist of the normal 18 Standard Model parameters and three new parameters: A, 
tan/5 = Vu/vd (the ratio of the vacuum expectation values of the neutral members of the Higgs 
doublet fields, and H^, responsible for giving mass to up-type and down-type quark fields, 
respectively) and sign(/i). This is the result after trading in the B and parameters for 
v"^ = vf^ + and tan/?, and using the observed value of mz to fix v. For this study, we have 
taken tan/5 = 2 and sign(yu) = —1. Further, we will consider only Nm = 1 models in this 
paper. 

The masses of the sparticles as a function of A, resulting from the above inputs, are easily 
summarized. For Nm = 1 we have the following approximate rules of thumb when A is of order 
100 TeV: 

mjo ~ 1.35GeVx A(TeV), mj+ ~ 2.7 x A( TeV) , m~ ~ 8.1 GeV x A( TeV) , 
m^^ ~ 1.7GeV X A(TeV), mj^ ~ 3.5 GeV x A( TeV) , m~ ~ 11 GeV x A( TeV) , (3) 

where Mi is the U(l)-gaugino soft-supersymmetry-breaking mass, m~ the gluino mass, m~ a 
typical u or d squark mass, and ^7^, ""^j^ are the masses of the right- and left-handed sleptons 
(all evaluated at the appropriate Q < ITeV). The masses of various sparticles in the case of 
Nm = 1 are plotted as a function of A (taking f = g = 1) in Fig. |I[ 

The range of A of interest in this model is easily determined using the dependence of the 
sparticle masses on A. An approximate upper limit is set by requiring that all superpartner 
masses be < 1.5 TeV. An approximate lower limit comes from requiring that the right-handed 
slepton be heavier than the rough lower limit of ~ 70 — 80 GeV from current LEP-2 data [§ 
(see later discussion). From Fig. |1] and Eq. (^, we see that A values below 50 — 60 TeV would 
yield too small a value for itlj^. For A values above about 150 TeV the squarks and gluino 
would have masses of order 1.5 TeV, which is somewhat beyond the range for which the model 
would provide a comfortable solution to the naturalness problem. 

Assuming that supersymmetry is spontaneously broken in the hidden sector [Q, Q , a mass- 
less Goldstone fermion, the goldstino, arises. Its coupling to a particle and its superpartner is 
fixed by the supersymmetric Goldberger-Trieman relation: 

A„t = -];/"5^G„ + h.c., (4) 

where j^" is the supercurrent, which depends bilinearly on all the fermion-boson superpartner 
pairs of the theory and Ga is the spin-1/2 goldstino field (with spinor index a). ^/F is the 
scale at which supersymmetry-breaking occurs in the hidden sector. When gravitational effects 
are included, the goldstino is "absorbed" by the gravitino (^3/2), the spin-3/2 partner of the 
graviton. By this super- Higgs mechanism, the goldstino is removed from the physical spectrum 
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Figure 1: m~o, m~+, m~ , m~ , m~, mr and m~ as functions of A in the Nm = 1 GMSB 
scenario, taking f = g = 1 and = I.IA. 

and the gravitino acquires a mass {1^3/2). In models where the gravitino mass is generated at 
tree-level one finds 

where Mp is the reduced Planck scale [Mp = {SttGn)'^^^ ~ 2.4 x lO^^GeV]. The helicity 
±1 components of the gravitino behave approximately like the goldstino, whose couplings to 
particles and their superpartners are determined by Eq. (^. In contrast, the helicity ±| com- 
ponents of the gravitino couple with gravitational strength to particles and their superpartners, 
and thus can be neglected. It is convenient to rewrite Eq. (13) for the gravitino mass as follows: 




In the models we consider, the gravitino will be the lightest supersymmetric partner. The 
Xi and Ti will be substantially heavier; whichever of these two particle is the lighter will be 
called the next-to-lightest supersymmetric particle (NLSP), and will decay to the gravitino and 
its SM partner. If R-parity is conserved, as we assume, the gravitino will then be stable and 
thus will be a candidate for dark matter. If the gravitino LSP is too heavy > few keV) its 

relic density will overdose the universe in most cosmological scenarios FromEq. d), this 
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means that v F values above roughly 3000 TeV are disfavored. On the other hand, there is no 
particular problem with very light gravitinos (eV masses). Although these will not contribute 
significantly to the total mass density of the universe, it may turn out that the main component 
of the dark matter has another source. Some examples are: the QCD axion, its supersymmetric 



partner (the axino 0) or the lightest stable particle in the GMSB messenger sector ||10|| . 

The scale of supersymmetry breaking, a/F, is a very crucial parameter for the phenomenol- 
ogy of GMSB models. However, it is highly model-dependent. In hidden-sector models, values 
of \/F ~ lO^-lO'' TeV are required in order to generate sufficiently large supersymmetry- 
breaking in the sector of MSSM fields [|TT|, [^]. In particular, one can derive [|rT| the approximate 
inequality 

' ^ '>f(A)(^), (T) 



2000 TeV / - ' VlOOGeV/ \27.5y/K 
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where / = F/Fs ~ ~ 2.5 x 10"^ /g^, with being the coupling of the gauge 

group responsible for communicating (via two-loop diagrams) supersymmetry breaking from 
the hidden sector to the messenger sector. Perturbativity would require (yf^ < 1. For mj^ > 

75 - 80GeV (the rough LEP-2 limit) and / = 2.5 x 10^ Eq. (0) yields > 10^ TeV 
{\/F > 5000 TeV) for = 1 (A^m = 4). Thus, even allowing for the roughly factor of 2 or 3 
uncertainty in the approximate lower bound, a/F must lie near its upper limit from cosmology. 
In direct-transmission models, the two-loop communication between the hidden and messenger 
sectors is eliminated, / in Eq. (|^ is effectively of order unity, and \/F can be as low as 100 TeV 
in phenomenologically viable models. 

In what follows, we will focus on the (A^ = 1) model in which the x? is the next-to- 
lightest supersymmetric particle. Q We will employ the specific input boundary conditions given 
earlier. With these inputs, the model is completely specified in terms of just two independent 
parameters, A and ^/F. As explained above, A determines the gaugino and sfermion masses, 
while \/F determines the properties of the gravitino and (see below) the lifetime of the Xi- 
For ^/F in the 100-3000 TeV range, the small size of the couplings, Eq. (^), imply that all the 
sparticles other than the Xi undergo chain decay down to the Xi- The Xi finally decays to a 
two-body final state containing the ^3/2: Xi ~^ 193/2 is the only allowed two-body mode for 
A < 80 TeV, and this mode remains dominant out to the highest A values considered. The 
CT for Xi 193/2 decay takes the form: 



,100 GeV , , V ^ , 

cr ~o z ~ -130 — um 8) 

V ^x;=B-793/2 1 mzo I 100 TeV ' ^ 



X 



^ Focus on this model was originally motivated by the 6677 event at the Tevatron p9[ |. However, recent 
Tevatron analyzes pc| ] rule out most, if not all, of the parameter space for which this model could have explained 
this event. 

^The other two-body modes that become kinematically allowed at high A are Xi ~^ ^93/2 and Xi — > h^9?,/2- 
However, for a Xi that is nearly pure bino, the x\ ~* ^93/2 width is suppressed relative to that for Xi ~^ 753/2 
by a factor of [s'm 9w / cos 9w]^ and the Xi ^ ^°53/2 partial width is proportional to the square of the very 
small higgsino component of the Xi- 
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If ^/F ~ 3000 TeV (the upper limit from cosmology), then cr ~ 100m for m~o = 100 GeV; 

•s/F ~ 100 TeV implies a short but vertexable decay length. Thus, the signatures for this 
GMSB model are crucially dependent on ^/F . In particular, for large / = F/Fs, and hence 
large the neutralino decay is most naturally characterized by a substantial decay length, 
possibly of order tens to hundreds of meters. Events in which one or more of the neutralinos 
travels part way through the detector, and then decays, can be a substantial fraction of the 
total, as can events in which all the neutralinos exit the detector before decaying. Thus, at the 
very least, it is highly relevant to assess how our ability to discover supersymmetry changes 
as a function of \/F. In the following sections of this paper we will focus on this issue in the 
context of the DO detector at the Tevatron collider. 

Before proceeding with the Tevatron analysis, let us return to the question of what the 
current LEP-2 limits are upon the simple GMSB model we are considering. The latest limits 
[^[including ^/s = 183 GeV running are quite significant, at least in the limits of either small 
yF (for which Xi ~^ 793/2 is prompt) or large \/F (for which the Xi is invisible). In the 
small \/F /prompt-photon limit, slepton mass limits of > 75 — 80 GeV arise from looking for 
^bIr events where Ir ixi ~^ ^193/2- Further, m~o > 80 GeV is required by the absence of 
excess 77 events. In combination, these require A > 60 TeV (see Fig. |T]). In the large 
a/F /invisible-Xi limit, the normal SUSY search for Enin with Ir — > £xi implies rrij^ > 80 GeV, 

implying A > 55 TeV. The case of intermediate ^/F has not yet been specifically analyzed, 
but we would anticipate that (conservatively) A < 50 — 55 TeV would be excluded by such an 
analysis. Thus, as stated earlier, we shall focus our analysis on the A > 50 — 60 TeV region. 

Finally, we wish to caution that many alternative GMSB models have now been con- 
structed. In particular, the various direct-transmission models are very different from one 
another and need not have the simple M^-scale boundary conditions of Eqs. (|I]) and (H). The 
relations between sparticle and gravitino masses vary tremendously. As an extreme example, 
in the model of Ref. the gluino is the LSP and the gravitino is sufficiently heavy as to be 
phenomeno logically irrelevant. 

The remainder of the paper is organized as follows. In Sec. 2, we describe in detail the phe- 
nomenology of the x?-LSP model and give results for all the signals considered. In particular, 
we describe in detail the delayed Xi decay signals and illustrate the regions of model param- 
eter space for which they would allow detection of a supersymmetry signal. Sec. 3 presents 
discussion and conclusions. 



2 Collider phenomenology 

In the Introduction, we emphasized that our ability to detect supersymmetry at the Tevatron 
will depend substantially on ^/F. The case where VF is sufficiently small that all Xi decays 
are prompt received early attention H^, |1^, ^, |15|, |16[. Of these references, only ||15[ performs 
a quantitative Tevatron collider study with the full boundary condition constraints of Eqs. (0) 
and (0). The results of Ref. |15[ were confirmed in Ref. [jl8[. In these latter analyzes, the 
prompt-decay hadron collider signal for supersymmetry considered was events containing two 
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or more isolated photons, deriving from decay of two or more neutralinos, plus missing energy 
(and possibly other particles). The case where ^JF is very large and most x? decays occur 
outside the detector was first analyzed in Ref. ITS; it is equivalent to conventional super- 



symmetry phenomenology (in which, for example, one looks for jets plus missing energy or 
three leptons plus missing energy) with the constraints implied by Eqs. (|I]) and (Q) among 
the sparticle masses. For moderate a/F, it is most appropriate to look for signals that would 
be characteristic of events in which the Xi decay is delayed, but still occurs within the de- 
tector (or within a region that could be probed by a simple addition to the detector). In 
Ref. [^] one signal of this type was examined in the context of the DO detector: namely, the 
'outer-hadronic-calorimeter' (OHC) signal that results when the delayed decay takes place in 
one of the outermost DO hadronic calorimeter cells. Here, we will consider in addition, two 
other delayed-decay signals. The first is the 'impact parameter' (6) signal. If the (transverse) 
impact parameter of a photon appearing in the electromagnetic calorimeter could be shown to 
be non-zero, this would constitute a signal for a delayed decay. The extent to which the reso- 
lution for such an impact parameter measurement by the CDF and DO detectors is adequate 
for this signal to be useful is an important issue. The second is the 'roof-array' (RA) signal. 
For large many x5 decays will occur outside the muon chambers. A two-layer detector 
placed on the roof of the DO detector building could pick up a portion of such delayed photons 
and possibly provide a dramatic signal. In fact, we will demonstrate that sensitivity to GMSB 
models with substantial Xi decay length will be maximized by implementing all three (OHC, 
6, RA) delayed-decay signals. 

We will now outline the results we have obtained for these signals and the associated 
procedures. Events were generated simultaneously for all SUSY production mechanisms by 
modifying ISASUGRA/ISAJET p4[ to incorporate the GMSB boundary conditions, and then 
forcing delayed x? decays according to the predicted cr. We used the toy calorimeter simulation 
package ISAPLT. We simulated calorimetry covering |?7| < 4 with a cell size given by Ai? = 
A?7 X A0 = 0.1 X 0.0875 and took the hadronic (electromagnetic) calorimeter resolution to 
be Q.l/\fE (O.IS/V^). The DO electromagnetic calorimeter was simplified to a thin cylinder 
with radius r = Im and length —2m < z < +2m. Also important to the analysis was the 
central outer hadronic DO calorimeter (OHC), which we approximated as occupying a hollow 
solid cylinder defined by —2m < z < +2m and radial region 2m < r < 2.5m. It is important 
that the DO OHC is segmented in Ai? and that there are also several layers of inner hadronic 
calorimeter. 

The signals are defined in terms of jets, isolated leptons, isolated prompt photons and 
isolated delayed-decay photons {i.e. emerging from x? decays occurring with substantial delay). 

• A jet is defined by requiring |?7jet| < 3.5 and > 25 GeV (for Ai?coai. = 0.5). 

• An isolated lepton is defined by requiring \ri^\ < 2.5 and > 20, 15, 10 GeV for the 
most energetic, 2nd most energetic and 3rd most energetic lepton, respectively. f\ The 
criterion for isolation is Et{AR < 0.3) < 4 GeV (summing over all other particles in the 
cone surrounding the lepton). 



^We never consider more than three isolated leptons, and all leptons are prompt. 
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• A photon is a prompt-photon candidate if it emerges from a Xi decay that occurs before 
the Xi has reached the electromagnetic calorimeter. The {rj^, 4>^) of such a photon is 
defined by the direction of the vector pointing from the interaction point to the point 
at which it hits the electromagnetic calorimeter. (This is generally not the same as 
the direction of the photon's momentum.) An isolated prompt photon is defined by 
requiring jr^^l < 1 and Ej< > 12 GeV, with isolation specified by Et{AR < 0.4) < 4 GeV, 
summing over all other particles in the cone surrounding the point at which the photon 
hits the EM calorimeter. Photons that emerge within the hollow cylinder defined by the 
electromagnetic calorimeter, but that are not isolated, are merged with hadronic jets as 
appropriate. 

• An isolated delayed-decay OHC signal photon is defined as one that emerges from a 
Xi — > 7^3/2 decay that takes place inside an OHC cell and that has E'J > 15 GeV, with 
isolation specified by Et{AR < 0.5) < 5 GeV, summing over all other particles in the 
cone surroimding the location of the Xi decay. This implies, in particular, that there 
should be very little energy deposited in the inner hadronic calorimeter cells in the same 
(?7, 0) location as the OHC cell in which the Xi decay occurs. 

• An isolated delayed-decay impact-parameter signal photon is defined as one that emerges 
from a. Xi ^ 193/2 decay (occurring before the Xi has exited the tracking region) and 
that passes through the central electromagnetic calorimeter (with the dimensions and 
physical location specified earlier). The photon is required to have Ej. > 15 GeV. No 
isolation requirement is imposed. Nor is there any specific requirement on the apparent 
rapidity {i.e. as measured by the location at which it passes through the electromagnetic 
calorimeter) . 

• A delayed-decay roof-array signal photon is defined as one that emerges from a Xi 
7^3/2 decay that occurs at a vertical distance in the range 6.5m < z < 16.5m from the 
interaction point and has Ej, > 15 GeV. This means that the Xi will have decayed past 
the muon detector system but below the roof of the DO detector building. The delayed- 
decay photon is required to pass through a 28m x 38m rectangle centered at 16.5m above 
the interaction point. This is the approximate size of roof array that could be positioned 
so that no part of the array would have a large amount of material between it and the 
6.5m height referred to earlier. No further isolation requirements or other cuts are placed 
on the photon. 

We now give the specific requirements for each type of signal. First, for all signals, events 
are retained only if at least one of several sets of reasonable trigger requirements are satisfied. 
The triggers considered were the following. 

1. I^T > 35 GeV. 

2. four or more jets with E^^?* > 15 GeV. 

3. (a) two or more jets with E^^ > 30 GeV; (b) > 40 GeV. 
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4. (a) one or more leptons with > 15 GeV, \ri£\ < 2.5 and isolation defined by Et{AR < 
0.3) < 5 GeV; (b) > 15 GeV. 

5. (a) two or more leptons with > 15 GeV, \T]i\ < 2.5 and isolation defined by Et{AR < 
0.3) < 5 GeV; (b) > 10 GeV. 

6. two or more photons with E^ > 12 GeV, |?7^| > 1 and isolation specified by Et{AR < 
0.4) < 4GeV. 

The additional signal-specific cuts are itemized below. 
• The standard jets-plus-missing-energy signal. 



We employ DO cuts ||2T|: (a) riQets) > 3 — labelled = 1, 2, 3 according to decreasing Et; 
(b) no isolated (^^^'^■(Ai? < 0.3) < 5 GeV) leptons with Et > 15 GeV; (c) > 75 GeV; (d) 



0.1 < A(f){]^Tjk) < vr - 0.1 and ^(A0(^t, Ji) - vr)2 + (A^d^T, ^2))^ > 0.5. The background 
cross section level has been estimated by DO to be ctb = 2.35 pb. The signal will be deemed 
observable if: (i) there are at least 5 signal events; (ii) (Ts/ctb > 0.2; and (iii) Ns/^/Nb > 5, 
where Ns and Nb are the numbers of signal and background events. 

• The tri-lepton signal. 



We employ the cuts of Ref. |2^: three isolated leptons, riQets) = 0, > 25 GeV and 
\M{iiij)—mz\ > 8 GeV (z 7^ j). For these cuts, the cross section for the sum of all backgrounds 
to the tri-lepton signal is 0"^ = 0.2 fb. 

• The tri-lepton-plus-prompt-photon signal. 

We impose the same cuts as for the tri-lepton signal, but require in addition that there be at 
least one prompt photon (as defined above) present. We shall presume that the background to 
this signal will be negligible. This seems reasonable given the already quite small background 
to the simple tri-lepton signal. 

• The two-prompt-photon-plus-missing-energy signal. 

Following Ref. [|l5l, we require: (a) at least two isolated prompt photons (as defined above); 



and (b) > W""- Detection efficiency of 80% (100%) is assumed if < 25 GeV (> 25 GeV). 
This signal should be completely free of background provided is adjusted appropriately 
as a function of luminosity. The signal will be deemed observable if there are 5 or more events. 

• The delayed-neutralino-decay signal resulting when the Xi decay takes place in the out- 
ermost hadronic calorimeter cell of the DO detector. 
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In an event in which a Xi 7^3/2 decay occurs inside one of the outer hadronic calorimeter 
(OHC) cells, the 7 will deposit all its energy in the cell. By demanding substantial 7 energy and 
isolation (precise requirements were given earlier) for this deposit, along with other criteria, 
backgrounds can be made small. More specifically, we require that the event fall into one of 
three classes defined by the following sets of requirements: 

1. (a) ra(jets) > 3; (b) at least one delayed-decay OHC photon; (c) J^t >-^™™. 

2. (a) any number of jets; (b) two or more delayed-decay OHC photons; (c) >-^™™. 

3. (a) any number of jets; (b) at least one prompt photon; (c) at least one delayed-decay 
OHC photon; (d) >^t'°- 

In the absence of the needed detector-specific study, we will assume zero background to 
1.+2.+3. for the same values that will be employed for the prompt-27 signal. Ob- 
servability of the delayed-decay OHC signal is assumed if the number of events for 1. +2.-1-3, 
is 5 or more. Of the above three sub-signals, 1. provides by far the largest event rates in the 
moderate \/F region. 

• The impact-parameter signal. 

We require riQets) > 3, > and one or more delayed-decay impact-parameter 

photons (as defined above) with transverse impact parameter b larger than 2cm. The DO 
detector in Run-II will have a pre-shower installed and is expected to achieve a la resolution 



of about 0.2cm in the transverse impact parameter. Thus, our 2cm requirement amounts to a 
lOcr deviation from zero impact parameter. 

• The roof-array signal. 

We require n(jets) > 3, > $t^^ and one or more delayed-decay roof-array photons. 
If there is directionality or timing for the roof array, additional cuts that ensure tiny back- 
grounds could be incorporated without significant losses in the signal event rate. However, 
such additional cuts are not imposed in this study. 

We will now summarize our results for all these channels. A limited subset of these results 



appeared in Ref. ||T^. In the present paper, our discussion will focus on the luminosities 
expected in Run-II and at TeV33. We begin with the jets-plus-missing-energy and tri-lepton 
signals. 

^Our additional criteria will be chosen so that the cross section for producing an isolated energetic long-lived 
kaon in association therewith is small. Further, any such kaon will interact strongly and is almost certain to 
be absorbed before reaching the OHC. 

^We will not discuss the delayed-decay signals for Run-I luminosity here. The OHC signals give only a 
marginal increase in the range of parameter space for which a supersymmetry signal might emerge for the low 
L ~ 100 pb~^ integrated luminosity accumulated in that run. Only if the stringent cuts we employ for the 
delayed-decay OHC signals could be weakened without encountering backgrounds would present data allow 
elimination of a significant additional portion of parameter space. Background studies are now in progress [2^ . 
The impact-parameter and roof-array signals cannot be implemented for the Run-I data set. 
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• The jets-plus-missing-energy-signal. 

Employing the cuts specified above, for L = 0.1 fb~^ (Run-I), we found in Ref. [|l^ that 
this signal was viable for any \/F if A < 25 TeV. For both L = 2 fb~^ and L = 30 fb~^, 
the signal is viable for any \/F if A < 30 TeV. Given the lower bound of ~ 55 TeV 
coming from the lower bound on fnj^, it would appear that the jets-plus-missing-energy 
signal will not be useful at the Tevatron. 

• The tri-lepton signal. 

Contours of constant cross section for this signal are presented in Fig. |^. Comparing to 
the background cross section, as = 0.2 fb, we see that in order to achieve S/\fB = 5, 
we require a signal cross section of a = 8.4 fb at L = 0.1 fb~^ (Run-I) a = 1.6 fb at 
L = 2 fb~^ (Run-II) and a signal cross section of cr = 0.41 fb at L = 30 fb^^ (Tev33), 
respectively. Fig. |^ then shows that the parameter space for which supersymmetry can 
be discovered using this signal extends to A ~ 50 TeV for L = 0.1 fb~^, to A ~ 65 TeV 
for L = 2 fb~^ and to A ~ 75 TeV for L = 30 fb~^. Thus, the tri-lepton discovery 
mode will allow supersymmetry detection out to substantially higher A than does the 
jets plus missing energy mode. (This is to be expected since the tri-lepton signal tends 
to be stronger than the jets-plus-missing-energy signal at Tevatron energies ||2^.) In 
particular, it can probe above the ~ 55 TeV lower bound on A coming from LEP-2 limits 
on mj . 

For both these signals, we should note that if a/F is small enough then, in association with 
the jets or leptons, one will detect prompt photons (coming from the decays of the Xi's that 
are inevitably present). The prompt-two-photon signal to be discussed below will also be 
in evidence. There will then be little doubt that nature has chosen a GMSB model rather 
than a SUGRA |^ model with similar boundary conditions for the gaugino and sfermion soft- 
supersymmetry-breaking masses. However, if a/F is large there will be no prompt photons 
and no hint for the GMSB character of the model from these two signals. To illustrate these 
points, we turn to the prompt-photon signals. 

• The tri-lepton-plus-prompt-photon signal. 

Contours of constant cross section for this signal are presented in Fig. ^. The solid 
contours are the same as plotted in Fig. The dotted contours correspond to 5 events 
at Tev33 (a = 0.16 fb at L = 30 fb"^) and 5 events in Run-II (a = 2.5 fb at L = 2 fb"^). 
Five events should be adequate for discovery given the negligible background expected. 
We observe that prompt photons will be seen in association with the tri-lepton signal 
only for < 600 TeV at L = 2 fb"^ and only for ^/F < 1600 TeV at L = 30 fb"\ 

We now turn to the impact-parameter (6), outer-hadronic-calorimeter (OHC), roof-array 
(RA) and prompt-two-photon (27) signals. In all these cases, we shall presume that our cuts 

^As usual, the SUGRA acronym refers to models in which supersymmetry breaking in a hidden sector is 
transmitted to the supersymmetric partners of the SM particles by gravitational interactions rather than gauge 
interactions. The gravitino in such models is heavy and an LSP x5 stable in the absence of R-parity violation. 
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Figure 2: We present cross section contours in fb units for the tri-lepton signal (as defined in 
tfie text). 

can be chosen so as to ehminate any background. In particular, we will consider cross sections 
and event rates for the cuts stated earlier with various choices of A complete summary 
of our simulations is contained in Figs. ^, ^ and ^ where we give cross section contours for 
the above signals in the (V^, A) parameter space for^™° = 30, 50 and 70 GeV, respectively. 
We now discuss the implications of the results contained in Figs. ^ ^ and ^. 

• The prompt-27 signal. 

Based on the background results of |^ (see also |^) and eye-ball extrapolations thereof. 



we estimate that our jet requirements and photon energy/isolation requirements are such 
that = 30, 50, 70 GeV is sufficient to eliminate all background for luminosities of 
L = 100 pb~^, 2 fb~^, 30 fb~^, respectively. Thus, to ascertain the discovery reach of this 
channel for L = 100 pb~^, one refers to the a = 50 fb contour of Fig. ^; for L = 2 fb~^, 
one refers to the a = 2.5 fb contour of Fig. ^; and, for L = 30 fb~ , one refers to the 
(T = 0.16 fb contour of Fig. ^. From these contours, we conclude that the signal is viable 
in a region described as follows. For L = 2 fb~^, the region extends from small a/F up 
to ~ 600 TeV if A < 35 TeV, increasing to ~ 1000 TeV at A ~ 80 - 100 TeV; 
the signal is not viable beyond A ~ 100 TeV. For L = 30 fb~^, the region extends 
from small \/F up to \/F ~ 900 TeV if A < 35 TeV, increasing to ^/F < 1000 TeV at 
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Figure 3: We present cross section contours in fb units for the tri-lepton-plus-prompt-photon 
signal (as defined in the text). In addition to the cuts for the tri-lepton signal, we impose the 
requirement that there be at least one prompt photon present. The solid contours are the same 
ones presented in Fig. ||. The dotted contours correspond to 5 events at Tev33 {L = 30 fb~^) 
and 5 events in Run-II (L = 2 fb~^). 

A ~ 120 - 130 TeV; the signal is not viable beyond A ~ 150 TeV. 
• The OHC signal(s). 

We provisionally adopt the same^^^'^ cuts of^f'' = 30, 50, 70 GeV for L = 100 pb"\ 
2 fb~^, 30 fb~^, respectively, as employed for the prompt-27 signal. Assuming that 5 
events are sufficient {i.e. that the background is negligible for our other strong cuts on 
jets and photon energy/isolation) to give an observable signal, we look at the same a 
contours as a function of L as in the prompt-27 case, with the following results. For 
L = 2 fb~^, the OHC signal provides a useful, but not enormous, extension of the region 
of parameter space for which supersymmetry can be detected as compared to the portion 
of parameter space covered by the prompt-two-photon signal or jets-plus- missing-energy 
signal. This extension is mainly in the \/F < 800 — 1000 TeV, A < 50 TeV region. Of 
course, we noted earlier that LEP-2 data probably already excludes this A region. For 
L = 30 fb~^, the reach of the OHC signal is greatly extended; supersymmetry detection 
may be possible via the OHC signal for essentially all of the \/F < 1500 — 2000 TeV, 
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Figure 4: We present cross section contours in the {■sfP, A) parameter space for the (a) impact- 
parameter (6), (b) outer-hadronic (OHC), (c) roof-array (RA), and (d) prompt-two-photon (27) 
signals. Contours are given at cr = 0.16, 2.5, and 50 fb (the outermost contour corresponding 
to the smallest cross section, etc.). We note that a = 0.16,2.5,50 fb corresponds to 5 events 
at L = 30,2,0.1 fb~^, respectively. Hence the labels Tev33, Run-II and Run-I, respectively. 
In this figure we have taken = 30 GeV for all signals. Other jet and photon requirements 
are as specified in the text. 
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Figure 5: As in Fig. g but for = 50 GeV. 
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Figure 6: As in Fig. g but for = 70 GeV. 
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A < 120 TeV portion of parameter space for which the prompt-two-photon signal is not 
viable. Conversely, supersymmetry detection at the \ow-\/F points at large A for which 
the OHC signal falls below 5 events would be possible via the prompt-27 channel. In 
addition, there is a substantial region in which both these signals are viable. Since the tri- 
lepton signal only extends out to about A < 75 TeV for L = 30 fb~\ it is apparent that, 
at the Tev33 luminosity, the prompt-27 and OHC signals combine to significantly extend 
to higher values of A the region where supersymmetry can be detected. In particular, 
A values substantially beyond the approximate limit of A > 55 TeV coming from LEP-2 
data can be probed. 

We make two further remarks regarding the OHC signals. 

— Of the three delayed-decay OHC signals, the first (> 3 jets -|- > 1 OHC 7) is 
observable (i.e. yields at least five events) for the entire (V^, A) region where the 
OHC signals are shown to be viable. The second (> 2 OHC 7's) yields > 5 events 
for only a very small portion of this region. The third (> 1 prompt 7 -|- > 1 OHC 
7) yields > 5 events for the lower (roughly, half) \/F portion of the region for which 
the combined OHC signals are viable. 

— If simulations eventually show that the first OHC signal is background-free for 
smaller^™™ than employed here (as we are hopeful will be the case), the portion 
of parameter space for which it is viable expands; e.g. using = 50 GeV at 
L = 30 fb~\ the > 5 event region would extend to ~ 1800 - 2300 TeV for 
A < 120 TeV. 

• The impact-parameter signal. 

Requiring^™" = 50, 70 GeV at L = 2,30 fb~^, respectively. Figs. || and | show that 
the b signal has an even larger region with > 5 events than does the OHC signal. For 
L = 2 fb~^, one finds 5 or more events for ^/F < 1200 — 1700 for A < 115 TeV. For 
L = 30 fb"\ one finds 5 or more events for VF < 2200 - 3000 TeV for A < 150 TeV. 
Only in a small corner with small ^/F and large A does the b signal fail, and this region 
is handily covered by the prompt-27 signal. As compared to the tri-lepton signal, the 
impact-parameter signal extends to substantially higher A values for \/F < 1600 TeV at 
L = 2 fb~^ and for y/F < 2700 TeV for L = 30 ih~\ 

From comparing Figs. ^ ^ and ^ we can see that weakening the L = 30 fb^^ missing 
energy cut to-^™™ = 50 GeV results in a slight increase in the > 5 event parameter space 
region. Weakening the cut still further to ^^'"^ = 30 GeV appears (we cannot be precise 
since our scan did not extend to high enough a/F) to result in a substantial increase to 
higher a/F of the > 5 event region. 

• The roof-array signal. 

For L = 2 fb-\ the RA signal has > 5 events up to /F < 1500 TeV for A < 40 TeV. 
In all of this region, the tri-lepton and impact-parameter signals, but not the tri-lepton- 
plus-prompt-photon signal, will also be viable. For L = 30 fb~ , the RA signal has > 5 
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events for large \fF < 2800 - 3000 TeV when A < 120 TeV. This A reach at high VF 
is substantially greater than achieved by the tri-lepton signal. As \/F decreases, the A 
extent of the > 5 event region decreases. 

From comparing Figs. |, ^ and | we can see that weakening the L = 30 fb"^ cut tol^^'"" = 
50 GeV, or, if possible without encountering backgrounds, = 30 TeV, substantially 
increases the high ^/F reach of the > 5 event region for A < 120 TeV. [| 

One especially important point regarding the delayed-decay signals discussed above is the 
following. If a/F is large, they are the only signals that will allow one to ascertain that 
nature has chosen a GMSB model, even if supersymmetry is detected in the jets-plus-missing- 
energy and/or tri-lepton modes. In particular, we saw above that the tri-lepton-plus-prompt- 
photon and prompt-two-photon signals will not be in evidence at large \/F, regardless of 
the magnitude of A. Without the delayed-decay signals, a SUGRA model with GMSB-like 
boundary conditions for the soft-supersymmetry-breaking masses of the gauginos and sfermions 
would then be indistinguishable from a GMSB model with a large value for a/F. 

The dependence of our signals on various energy variables may be useful in designing 
triggers and considering backgrounds. To this end, we present a figure showing the distributions 
in Ej, E^r^^ and^T for the impact-parameter and roof-array signals. These distributions are 
presented for the parameter space point a/F = 1400 TeV and A = 50 TeV; they are fairly 
representative of the results obtained at other parameter space points. In these distributions, 
we have imposed a more limited set of cuts. 

• In the case of the distributions presented for the impact-parameter signal, we require only 
that there be one or more delayed-decay impact-parameter photon (as defined earlier) 
with b > 2 cm, except that when plotting the E^ distribution, we remove the cut on E^. 

• In the case of the distributions presented for the roof-array signal, we require only that 
there be one or more delayed-decay roof-array photon (as defined earlier). 

In neither case do we impose jet cuts, except that when plotting da /dEj^, the jets must satisfy 
the definition of a jet as given earlier. In addition, no cuts are placed on If^T- In the case of 
the F;p* and E^ distributions, all jets and delayed-decay (impact-parameter in the first case, 
and roof-array in the second case) photons satisfying the minimum requirements noted above 
are included in the histograms. 

From the plots of Fig. |^ we observe the following. 

• Our E^ cuts accept essentially all events and could be strengthened somewhat without 
significant loss of signal event rate. 

• The F^?* spectrum is falling very rapidly, and it would be undesirable to strengthen the 
jet cuts unless it is absolutely necessary. 

^Since our parameter scan was limited to ^/F < 2800, we are miable to be quantitative in this statement. 
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Figure 7: For the impact-parameter and roof-array signals, we plot da/dE^, da/dE^^ and 

da/dpT for the parameter space point VF = 1400 TeV and A = 50TeV. The hmited cuts 
applied are described in the text. 
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Cuts significantly stronger than > 70 GeV will start to cause a significant loss of 
signal rate. Thus, we hope that such a stronger cut will not be needed to eliminate 
backgrounds. 




Figure 8: For > 50 GeV, we plot the same cross section contours for the impact-parameter 
signal as given in Fig. ^ (for b > 2cm), for the cases (a) b > 50cm and (b) b > 10cm. All other 
cuts are identical. 

It is also potentially important to know how sensitive the impact-parameter signal is to the 
impact-parameter resolution. To this end, we present in Fig. H the same cross section contours 
for the impact-parameter signal as given in Fig. ^ (for b > 2cm), for the cases (a) b > 50cm 
and (b) b > 10cm, keeping all other cuts (including = 50 GeV) the same. We observe 
that the b cut can be increased from 2cm to 10cm without significant damage to the signal. 
But, if the transverse impact parameter resolution were as poor as the 5cm-10cm resolution 
applicable during Run-I (for which the DO detector did not have a pre-shower device) a 
cut of 6 > 50cm would probably be needed to eliminate backgrounds and the region of viability 
of the impact-parameter signal would decrease very substantially. In fact, for all cases [{L = 
100 ph-\:^f'' = 30 GeV), (L = 2 fb-\^f?^'^ = 50 GeV) and (L = 30 {h-\:^f'' = 70 GeV)] the 
region for which > 5 events are found for the impact-parameter signal after requiring b > 50cm 
is quite similar to that for which > 5 OHC signal events are found. In contrast, the region for 
which > 5 impact-parameter signal events are expected for the b > 10cm or 6 > 2cm cuts is 
always substantially larger than the region for which > 5 OHC signal events are predicted in 
all the above three (L,^™™) cases. 
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3 Discussion and Conclusions 



We have studied a simple gauge-mediated supersymmetry breaking model in which the lightest 
neutralino, the Xi, is the next-to-lightest supersymmetric particle and decays to photon plus 
gravitino. The model is characterized by the parameters A (which sets the scale of the masses 
of the superpartners of the Standard Model particles) and ^/F (which specifies the scale at 
which supersymmetry is broken in the hidden sector and, thereby, determines the mass of 
the gravitino and the lifetime of the Xi)- We have argued that it is very possible (indeed, 
preferred in the simplest existing GMSB models) for the parameter to be sufficiently big 
that there is a large probability for the x? to decay a substantial distance from the interaction 
point. We have shown that the portion of the (\/F, A) parameter space for which a signal for 
supersymmetry can be seen at the Tevatron is greatly expanded by employing signals sensitive 
to such delayed decays. 

In particular, it is useful to compare (a) the standard SUSY signatures (jets plus missing- 
energy or three leptons plus missing-energy), (b) the prompt-two-photon signal and (c) the 
delayed-decay (impact-parameter, roof-array and outer-hadronic-calorimeter) signals. 

• Of the standard SUSY signatures, the tri-lepton signal is the stronger. Using it, discovery 
of GMSB supersymmetry will be possible for small to moderate A values, < 65 (< 
75TeV) at Run-II (Tev33), regardless of the size of 

• The prompt-two-photon signal is viable for small to moderate v^F, but extending to 
reasonably large A, A < 100 TeV (< 150 TeV) at Run-II (Tev33). 

• The delayed-decay signals are potentially viable in the large- a/F, large- A portion of 
parameter space not covered by the standard SUSY and prompt-two-photon signals. 
For Run-II, the delayed decay signals could allow supersymmetry discovery for all of 

< 1700 TeV, A < 100 TeV, extending at Tev33 to cover essentially the entire \fF < 
3000 TeV and A < 120 TeV region that is preferred in the context of the model we have 
explored. 

If \/F is large, the delayed-decay signals will be important even if A is in the small to 
moderate range (A < 65 — 75 TeV) for which the standard jets-plus-missing-energy and/or tri- 
lepton SUSY signatures will be visible. This is because the prompt-photon signals (tri-lepton- 
plus-prompt-photon or prompt-two-photon) will not be in evidence if a/F is big. Without 
employing the impact-parameter, roof-array and/or outer-hadronic-calorimeter delayed-decay 
signals there will be no hint that nature has chosen a GMSB model, as opposed to a SUGRA 
model with GMSB-like boundary conditions for gaugino and sfermion soft-supersymmetry- 
breaking masses. 

If delayed decay signatures are seen, the next goal will be to determine the parameters 
of the GMSB model. These include the all important value of \/F and the masses of the 

^In a more general GMSB model, it is certainly possible that the entire interesting region of parameter space 
might not be accessible using the delayed-decay signals. Still, we expect that the delayed-decay signals will 
always significantly expand the amount of GMSB parameter space at larger \/F for which a supersymmetry 
signal can be discovered. 
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supersymmetric particles. The latter will allow us to determine the overall scale A and the 
relative sizes of the soft-supersymmetry-breaking masses as a function of A. Important input 
information will include: (a) the distribution of x? decays as a function of distance from the 
interaction point, as reflected in the relative rates for different signals {e.g. the prompt-two- 
photon rate as compared to the roof- array rate); (b) the energy distribution of the photons 
from the Xi decays; (c) the time at which the photons from the Xi decays arrive at various 
elements of the detector {e.g. the electromagnetic calorimeter pre-shower, the outer-hadronic- 
calorimeter cell or the roof-array); (d) the missing-energy distribution; and (e) the types 
of different particles that appear in association with the supersymmetric events, and their 
energy distributions. All of these depend in a complicated way on the parameters of interest. 
For example, the distance distribution of Xi decays, the photon energies, and the timing of 
the photons from the decays all depend strongly on (cr)~o and m~o, and might allow their 
determination. f\ Even a rough measurement of these latter two quantities will allow us to 
estimate a/F using Eq. (||). Further, the energy distributions of the visible particles (jets, 
leptons, and photons) will reflect both the relative and absolute masses of the sparticles. For 
example, thresholds reflecting mj^ — m~o and/or mj^ — m~o might show up in lepton energy 
distributions. And so forth. Obviously, substantial statistics will be required for a relatively 
model-independent determination of \/F and the soft-supersymmetry-breaking masses. 

We note that all the delayed-decay signals we considered here are potentially useful for any 
model in which the decay of the Xi is not prompt. For example, in R-parity violating (RPV) 
SUGRA models the RPV couplings may be sufficiently small that the cr for Xi decay is quite 
substantial. The products emerging from the Xi decay (whether three leptons, one lepton 
plus two jets, or three jets) can emerge in the outer hadronic calorimeter, with large impact 
parameter, or outside the detector and pass through the roof array. There is also the possibility 
that R-parity violation could be present in GMSB models. For an appropriate balance between 
the magnitude of \/F and the size of the RPV couplings, the x? could have prominent decays 
both to RPV channels and to 7^3/2, and both types of Xi decays could be substantially delayed. 
In general, observation of the delayed-decay signals would be quite critical to determining the 
relative branching ratios for RPV and 7^3/2 decays of the Xi- A similar set of remarks would 
apply to RPV (SUGRA or GMSB) models in which a right-handed slepton is the lightest of 
the superpartners of the SM particles. Thus, the utility of the types of delayed-decay signals 
we propose will extend far beyond the NLSP=Xi, R-parity-conserving GMSB model context 
analyzed in detail here. 

Our study has focused on using the DO detector in the configuration planned for Run-II, 
with the addition of a (relatively cheap) roof-array detector. There are, of course, much more 

^As an example, one could measure the time of flight of the Xi to the OHC cell (by the time of the OHC 
energy deposit) and the energy of the deposit {i.e. the energy of the photon). The flight time is a function 
of cr, m~o and the average Xi energy. The photon energy is a function of the latter two. The correlation 
between the flight time and photon energy (as measured in a large number of events) might allow a roughly 
model- independent determination of ct and m~Q. 

^°The main difference, as compared to the Xi case, is that additional information would arise from observation 
of the visible track associated with a semi-stable charged Ir. Also, a long-lived might be absorbed inside 
the detector. 
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ambitious possibilities in which a decay volume and special detector, specifically designed to 
pick up delayed Xi decays, is constructed some distance from the main detector. Nonetheless, 
since the delayed-decay photon signals discussed here show great promise and since delayed 
decays in general could easily prove to be the hallmark signature of gauge-mediated supersym- 
metry breaking, not to mention R-parity violation, the CDF and DO detector groups should 
work to clarify and refine these signals and, in particular, the backgrounds thereto. 
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